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The phenomenology of heavy quarkonia production in hadron collisions is reviewed. 
The theoretical predictions are compared to data. Commonly used production models 
are shown to fail in explaining all the experimental findings. The shortcomings of these 
models are analysed and possible improvements are discussed. 
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1 Introduction 

The production of heavy quarkonium states in high energy processes probes the very border 
between perturbative and non perturbative domains. It provides therefore quite stringent tests 
of our understanding of Quantum Chromodynamics (QCD) and has hence recently attracted 
much of the theoretical and experimental interest. 

A large amount of data is now available, both from fixed target and from collider || [| 
experiments. In this talk I'll restrict myself to the inelastic hadroproduction case, and in partic- 
ular I'll focus on the collider results. This because that's where the most serious discrepancies 
between theory and experiments have shown up, demanding for a careful reanalysis of the 
quarkonium production mechanisms. 

The outline of the talk is as follows. I will first describe the general framework which is 
believed to encompass quarkonium production, and the way this framework has in the past 
been formalized into production models. The approximations which had been made will be 
pointed out, and the results of theoretical predictions based on these models will be compared 
to experimental data and found to be sometimes in disagreement. Possible reasons for the 
discrepancies will then be considered, and a new model which may help reconciling theory with 
experiments will be briefly presented and discussed. 

2 Production models 

When considering any bound state object production two key questions have to be answered: 

1. What is it made of? 

2. How are its components produced? 

As for the first question, there is evidence that heavy quarkonia are (mainly) composed of a 
heavy quark-antiquark pair, i.e. cc for charmonium (the J/tp family) and bb for the T family. 

As for the second questions, there seems to be general agreement that QCD is the theory of 
strong interactions. It is therefore within its context that the production of the heavy quark 
and antiquark has to be described. 

These two questions and the answers given above constitute the general framework to quarko- 
nium phenomenology. At this points the different models we can think of do diverge and provide 
different descriptions and hence different results. I'll now describe two models which have been 
put forward some years ago and widely used till today: Fritzsch's Color Evaporation Model [|J 
(CEM) and Berger and Jones' Color Singlet Model § (CSM). 

2.1 The Color Evaporation Model 

This model makes use of the parton-hadron duality hypothesis to relate the charmonium cross 
section to the quark-antiquark production cross section. In particular, it assumes that all 
the non perturbative effects that lead to the bound state formation cancel when considering 
inclusive final states, and therefore writes the production cross section for the quarkonium state 
H as: 




(1) 



On the right hand side we have the QCD cross section for producing a QQ pair, <j[QQ], 
integrated over the invariant mass range up to the mass tum of the lowest lying heavy-light 
meson (for instance the D meson in the charm case). The factor P# has to be determined 
phenomenologically and provides the only differentiation among the various quarkonium states. 
While being therefore quantitatively not very predictive, this model has its central feature in 
that the differential distributions (like energy, xf, or px dependencies) should be equal for all 
the quarkonium states considered. It is thus mainly on this ground that its degree of validity 
will have to be assessed. For a review of this model and its achievements see for instance ref . || . 



This model tackles the problem of quarkonium production from an opposite point of view with 
respect to the Color Evaporation Model. It tries to reach the highest possible predictivity at 
the price of making stronger approximations. 

Its central assumption is that the QQ pair that will form the bound state is produced by 
the hard (short distance) QCD interaction with the correct spin and angular quantum numbers 
for the quarkonium considered, in a color singlet (and therefore observable) configuration and 
with zero or small relative momentum: 



The non perturbative effects are assumed to factorize into a single parameter P n ^, which can 
be related to the wave function of the bound state (or its derivative, for P states) evaluated at 
the origin. It can be calculated within potential models (see for example ref. J/J) or extracted 
from experimental data of quarkonium decay widths, which can be calculated within the same 



The difference with the previous model can easily be appreciated: the structure of the 
quarkonium ("What is it made of?") is now defined much more precisely, and this allows 
detailed quantitative calculations for each different state H = n 2S+1 Lj. On the other hand, to 
assume that a quarkonium is well described by a QQ pair in a color singlet configuration may 
be too a crude approximation. 

3 Phenomenology 

When we make use of the CSM to obtain phenomenological predictions for quarkonia production 
we have to live with a double approximation. 

The first one is intrinsic in the definition of a quarkonium as a color singlet QQ pair with 
zero relative momentum: this amounts to neglect relativistic corrections, i.e. to work at lowest 
order in v 2 , v being the heavy quark speed in the quarkonium center of mass frame. 

The second is the perturbative order at which the cross section for producing the heavy 
quark- ant iquark pair is evaluated. The lowest order in a s was usually considered to represent 
the dominant term in the expansion. 

When calculating large pt production rates for quarkonia at the Tevatron, i.e. in pp collisions 
at \fs = 1800 GeV, within the above approximations we immediately find huge discrepancies 
between theory and experiment. The J/ip and ip' rates are found to be underestimated by one 
or two orders of magnitudes. 



2.2 The Color Singlet Model 



a[n 2S+1 Lj} = P nL a[QQ(n 2S+1 Lj,l)} 



(2) 




Figure 1: Theoretical prediction for the J ftp (a) and ip' (b) cross sections 
compared to CDF experimental data. 



We may at this point ask ourselves if are we operating within a wrong model or if we 
are just neglecting important contributions when making the aforementioned approximations. 
I'll pursue the second point of view and try to present the CSM as being the lowest order 
approximation of a wider model. It must therefore be supplemented by higher order terms, 
both in a s and in v 2 , in order to produce reliable predictions. 

Following the historical development I start considering higher order terms in a s . In 1993 
Braaten and Yuan || pointed out that the process 



" " \ 9 in (3) 
>-> J/ip gg 

while being of higher order (af vs. af) with respect to the lowest order one, 

gg -> J/tfj g (4) 

becomes however dominant at large px (say px > 5-6 GeV). The technical reason for this is 
that the process (0) is not suppressed by a form factor (m/px) 2 , m being the heavy quark mass, 
which is instead present in the cross section of The physical reason for the existence of this 
form factor is that it is difficult to produce the heavy quark and antiquark in a small spatial 
region, of size 1 /pt, and still keep their relative momentum close to zero so that they can form 
a bound state. This shortcoming is absent in the process @, as the gluon that fragments to 
the J/if) can have invariant mass of order m <C px and thus make the heavy quark- antiquark 
formation a longer distance process. Moreover, they could also show that in the limit of the 
gluon energy being much larger than its invariant mass the cross section could be approximated 
by the convolution of a kernel cross section for producing on-shell gluons with a fragmentation 
function for the gluon going to the J/i/j. 

The fragmentation functions for a gluon or a heavy quark to go into any quarkonium state 
can be calculated [§, ||] in perturbative QCD (with the exception of the usual non-perturbative 
parameters related to the bound state formation) and the cross section is given by 

d<r[H(p T )} _ f da[i(pT/z),jx\_ D jj . 



dpx J dpx 

In this equation o is the kernel cross section for producing the parton i, and Df is the frag- 
mentation function of i to the quarkonium H . fi is the factorization scale, to be taken of the 
order of px- 



Phenomenological predictions based on this formula [IU| have promptly shown that the CDF 
data on J/ip production could now be explained within a factor of two to five, depending 
on the input parameters, (see fig. |l|a) while predictions for ip' productions where still falling a 
factor of 30 below the data (see fig. |l|b). 

The solution to this discrepancy, which has been named the "CDF Anomaly", may have 
been found within the context I'm going to describe in the next Section. 



4 The Factorization Approach 



A rigorous frame for treating quarkonium production and decays has been recently developed 
by Bodwin, Braaten and Lepage [JTTj] . Their so-called "factorization model" expresses the cross 
section for quarkonium production as a sum of terms each of which contains a short- distance 
perturbative factor and a long-distance non perturbative matrix element: 

^=E^(0|Q?(A)|0> (6) 

F n are short-distance coefficients which can be calculated in perturbative QCD by expanding 
in powers of a s . A is a scale which separates short and long distance effects. The cross section 
is however independent of A as its effect is compensated by the A-dependence of the non- 
perturbative matrix elements. 5 n are related to the dimension of the operator O^. Finally, the 
matrix elements (0|O^(A)|0) can be defined rigorously in Non Relativistic QCD. They absorb 
the non perturbative features of the process and can either be extracted from data or calculated 
on a lattice. 

The main feature of this model, and the main difference with respect to the CSM, is that 
- pretty much like the Color Evaporation Model but in a much more sophisticated way - it 
takes into account the full Fock space structure of the quarkonium state. The latter is no more 
assumed to be represented by a color singlet QQ pair with the correct quantum numbers, but 
rather by an infinite series of terms: 

\H = n 2S+1 Lj) = 0(l)\QQ(n 2S+1 Lj),l) 

+ 0{v)\QQ{n 2S+ \L± l)j,,%)g) 

+ 0(v 2 )\QQ(n 2S+1 Lj,8)gg) + ... 

+ (7) 

The CSM can be recovered by taking the lowest order term in eq. 

Higher order components are suppressed by powers of v, but can become important if their 



associated short-distance coefficient F n in eq. (0) is large. Braaten and Fleming []T2| have 
shown that ip production via gluon fragmentation through a color octet 3 Si pair, being its 
short distance coefficient of order a s only, can easily overwhelm ordinary gluon fragmentation 
to a color singlet pair (of order a 3 ) and may thus explain the ip' abundance that CDF observes. 



5 Conclusions 

In this talk I have presented evidence that lowest order Color Singlet Model results are incapable 
of explaining CDF data of quarkonium production. Higher order terms both in a s and in v 2 are 
important and must be included to ensure a consistent and satisfactory description of quarkonia 
production. 

The interplay between higher orders in a s and in v 2 is far from being trivial, and quarkonia 
production has to be regarded as a two-parameter problem. A consistent framework for treating 
the double expansion can be found in the factorization model developed by Bodwin, Braaten 
and Lepage. 
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